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I. INTRODUCTION 



Although Quantum Chromodynamics (QCD), as the ultimate theory of strong interactions, should be able to give 
precise predictions of all physical quantities such as masses, mass splittings etc., its predictive power is obscured by 
technical difficulties. One way out is to employ effective models, which share some important features with QCD. 
The chiral quark model pj] (or semibosonized Nambu-Jona-Lasinio model) exhibits chiral symmetry breaking, which, 
together with confinement, is relevant for the low energy regime of QCD. Moreover it can be "derived" from QCD 
within the instanton liquid model of the QCD vacuum pi. 

The NJL model || has been therefore thoroughly examined against the low energy hadronic data and it turned 
out to give surprisingly good predictions for mass splittings [^J^] , electromagnetic properties || and axial properties 
of hyperons . The hadronic part of the isospin breaking effects which is due to Am = mj — m u mass difference is 
perhaps the most striking example of the accuracy reached by the model. In Ref. || the isospin mass differences have 
been calculated for the octet and decuplet of baryons. Although the error bars on these splittings are large (since the 
electromagnetic part of the splittings has to be subtracted) the model predictions fall nicely within the error bars for 
all spin 1/2 and 3/2 baryons for Am « 3.5 MeV. 

These results have been obtained within the perturbative approach, where one expands the collective hamiltonian 
describing baryonic states in terms of m s and Am up to terms linear in both parameters. It has been claimed in the 
literature pp| , pl| that, at least within the Skyrme model approach, the next orders in m s spoil the nice pattern of 
the isospin splittings obtained in the linear approximation. It is the purpose of this work to examine if this is also 
the case in the present model. As will be seen in the following, the inclusion of the terms of the order Am m s does 
not spoil the pattern of the isospin breaking, however Am which reproduces the experimental data is shifted towards 
the higher value of the order of 4.5 MeV. We consider this result as fully satisfactory. Other chiral models either 
underestimate the isospin splittings by factor of 2 or overestimate it by approximately the same amount J^] 

The crucial point of our analysis is the fact that though we are considering isospin breaking effects of non-strange 
current masses we are using the SU(3) version of the NJL model. In this framework the symmetry breaking operator 
is proportional to A 3 . In the SU(2) model it reduces to r 3 and most of the matrix elements after the collective 
quantization vanish. In SU(3) the polarization of the strange Dirac sea provides a number of terms which are crucial 
for the splitting pattern. 

II. MESONIC SECTOR 

The prominent feature of the NJL model consists in the fact that the same effective action describes meson physics 
and also, through the solitonic solutions, baryon physics. It is customary to fix the parameters of the model by fitting 
meson properties. This procedure leaves usually one free parameter, namely the constituent quark mass M. The 
regularization cutoff (or cutoff function) becomes then a function of M. Then M is fixed to fit one splitting in the 
soliton sector, eg. N— A, and then all other baryonic observables come out as predictions. In practice there is always 
some freedom in tuning other parameters like m s or Am, although they are always fairly constrained by the meson 
sector. 

There are in fact different ways to treat the meson sector. One can either solve the gap equations and then evaluate 
meson masses at zero momentum q — from the curvature of the effective potential. One can also evaluate the 
curvature at q 2 = — m 2 llcson . Or one can solve Bethe-Salpether equations for the meson propagators. Another way 
is to fix the parameters in the gradient expansion, which results in an effective meson lagrangians whose coupling 
constants are known from the meson scattering. Each of these methods produces slightly different results, which are 
not important for the gross features of the soliton sector, but may influence the numerics. This influence is not very 
significant if one is interested in the quantities which are not too small, however they might turn out important, if 
one considers such tiny effects as the isospin splittings. 

It is not our purpose to make a complete calculation of the isospin splittings. That would require to include the 
electromagnetic effects, which is certainly beyond the scope of this paper. Our goal is by far more modest: we want to 
investigate the effect of the strange quark mass on the hadronic part of the isospin splittings, very much in the spirit 
of the calculations performed in the Skyrme model |L^ |. Electromagnetic contributions will be simply parametrized 
by means of some convenient Ansatz; in the meson sector we use the Dashen Ansatz: 



1 See Introduction in Ref. [[Tl| 
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TO mcson = ™H + Q^- (!) 

where Q is the meson charge, c is a constant which does not depend on the meson in question and subscript H stands 
for the hadronic (quark mass dependent) part of the meson mass. 

We fix the model parameters in the meson sector in the first order in the quark mass matrix. In this way we avoid 
unnecessary complications due to the rf — rf — 7r 3 mixing. We also refrain from such problems as boson loops and the 
validity of the Dashen Ansatz. All these complications result in the uncertainty in the value of Am extracted from 
the mesonic data. 

The first order meson mass formulae with the Dashen Ansatz have been analyzed in Ref. || They imply: 

U = Am = 0.28. (2) 
m u + m d 

For m u + md = 12.2 MeV, which is required by the regularization prescription we use and which falls within the 
commonly accepted range of values for m u + md, we get Am = 3.4 MeV. Ignoring completely electromagnetic 
contributions would yield 1Z = 0.21 and Am = 2.6 MeV. 

It should be however stressed that full next-to-leading order analysis of mesonic masses within the chiral perturbation 
theory does not constrain Am at all Jll| . This is due to the lack of knowledge of the electromagnetic contribution 
to meson masses at next order in quark masses. It is however not excluded that at this order the Dashen Ansatz is 
violated. Recent analysis of the decay r; — > 37r suggests that 

(m|± - mg-o - ml± + m^) EM = 1300 ± 400MeV 2 , (3) 

in contrast to for the Dashen Ansatz. If this were true 1Z = 0.34 and Am = 4.2 MeV. 

Therefore we fix cutoff vs. M dependence as in Ref. ||, where the details can found. The experimental numbers 
which are used to fix the model parameters are f v = 93 MeV, m^ and m-K. Then it comes out that m u + md = 
12.2 MeV, m s ~ 150 MeV. The kaon decay constant is then constrained to /k = 105 MeV about 10% below the 
experimental value. The allowed ra ge for Am is then, as discussed above, 2.6 - 4.2 MeV, however the lower values 
are rather unlikely, since they come out by completely ignoring the electromagnetic contributions. 



III. SOLITONS AND THE QUANTIZATION OF ZERO MODES 

The effective Euclidean action for the semibosonized SU(3) NJL model reads (after integrating out the quark fields) 

S cff = -Sp ]n(-ip + m+ MU 15 ). (4) 

Here M is the constituent quark mass, m the current quark mass matrix and £/ 75 describes the 8 Goldstone modes of 
the SU(3) chiral symmetry. The real part of Eq.(|J) is assumed to be regularized with the proper-time regularization 
function of Ref. M and Sp denotes the functional trace. 

First we make use of the trivial embedding of G3 of the SU(2) chiral field Uq{x) — (ct( 2 ) + fysTTr)/ ' into the 
isospin subgroup of SU(3) according to 

£/(*)=( u ° ;). (5) 

The soliton solutions of SU(2) are also solutions for SU(3) and the embedding (|J) gives the correct constraint on the 
Hilbert space of the baryonic states |14| . 

The soliton solutions are found by employing the hedgehog Ansatz for the field Uq. The details are widely described 
in the literature. Following the treatment of Ref. p5| we quantize the soliton by introducing time dependent rotations 
of the hedgehog field: U(x,t) — A(t) U(x) A(ty. This rotation can be undone by rotating the quark fields: q — A(t)q 
and q = qA(tY . Defining "angular velocities" O: 

A^A = in E = l -\ a Q.% (6) 

one can rewrite S e s as: 

S cff = -Sp ln(<9 T + H + in E -i-y^mA). (7) 
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The following relation between Euclidean and Minkowski velocities holds: HIe — and Q E — Qe- 

Expanding (Q) up to the quadratic order in f2 (in Minkowski metric and in the chiral limit) one gets Ref. : 

L = -Ma + ^tt a I ab n b (8) 

where tensor of inertia I ab = diag(ii, I\, 1%, I 2 , h, I2, 12, 0) can be found in Ref. 0. 

To calculate the mass splittings one has to expand the effective action in powers of the current quark mass m = 
Mo ^0 - Ms ^8 - M3 ^3 (Ao = \/y3 1) with: 

Ho = —7=(m u + m<i + m s ), |ig = -i=(2m s -m u -m d ), ft = ^Am (9) 
Vo vl2 ^ 

where Am = ma — m u . 

Expanding the effective action up to terms of the order of m s , Am, m%, Am m s , m s fl and Am fl one gets: 
L m = -o m s + a ^m s D^ + ^ Amflfg'j , 

L m u = ~—^ m * D y 8 JK ab n b - AmDl a >K ab tt b , 
L m2 = Iml f JV (1 - D<§f + 3N ab D^D^ 
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Ar 



(jVbpSM - D$) + WaDgDW) , (10) 



where the constant <r is related to the nucleon sigma term S = 3/2 (m u + m d ) <r and = 1/2 Tr(At A Q AA 6 ). The 
mass spectrum obtained with the help of Lq + L m + L m Q was discussed in Refs. [QJ^]; there one can also find explicit 
formulae for K ab = diag(ifi, K±, K\, K 2 , K 2 , K 2 , K 2 , 0). Let us here only remind that the anomalous moments of 
inertia Ki are nearly entirely given by the valence part, whereas the contribution of the valence level to Ii amounts 
to approximately 60%. The quantities N ab — diag(iVi, Ni, Ni, N 2 , N 2 , N 2 , N 2 , Nq/3) have been derived in Ref. [||. 
Their values together with the values of I±, 2 and K\ t2 for different constituent masses are listed in Tab. | 

The lagrangian of Eq.([To|) reminds the one of the Skyrmion. The quantization proceeds as in the Skyrme model 
and has been described in detail in the literature Jlfj]. Let us here remind that at first one defines the quantities: 

J a — lab^b — HiD ib K ba ~ 6 a8 -^-j= (11) 

(i = 3 and 8, a, b = 1 . . . 8) which, as a result of the quantization procedure, are promoted to the spin operators J a . 
Note that the relation (O) depends on the quark masses. The wave function of the baryon state B — Y,T,Tq, J, J3 
belonging to the SU(3) representation TZ reads (see Appendix A of Ref. S): 



TZ, B > = VdimTZ (Y, I, h I (A) \ -Y' , J,-J 3 ) , (12) 



where the right hypercharge Y' is in fact constrained to be —1. The lowest SU(3) representations which contain states 
with Y — 1 are: TZ =8 and TZ —10. The quantized hamiltonian from Eq. (JsJ) reads: 

= M c i + H S \j( 2 ) + #su(3)) (13) 



su( 2) = — 0>(SU(2)), 



SU(3) 



1 

2^2 



C 2 (SU(3))-C 2 (SU(2))-^ 



Here C 2 denote the Casimir operators of the spin SU(2) and flavor SU(3). M c i is the classical soliton mass. It has 
been calculated by many authors and its value turns out to be relatively large: M c i ~ 1.2 GeV. This is a common 
problem for all chiral models. There are however some negative corrections to it, like Casimir energy or rotational 
band corrections which might bring M c i to the right value. In this paper, instead on insisting on the calculation of 
the absolute masses, we will concentrate on the mass splittings. 
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TABLE I. Moments of inertia for different constituent masses 



M 


E [SU(2)] 


h 


h 


Ki 


K 2 


N 


Ni 


N 2 


[MeV] 


[MeV] 


H 


[fm] 


[fm] 


[fm] 


[fm] 


[fm] 


[fm] 


363. 


60.32 


1.512 


0.720 


0.606 


0.372 


0.765 


0.647 


0.496 


395. 


58.14 


1.285 


0.618 


0.438 


0.290 


0.704 


0.500 


0.408 


419. 


56.14 


1.178 


0.569 


0.369 


0.255 


0.668 


0.438 


0.370 


423. 


55.52 


1.156 


0.560 


0.357 


0.250 


0.658 


0.426 


0.362 


442. 


52.52 


1.070 


0.521 


0.315 


0.229 


0.603 


0.379 


0.329 



The hamiltonian up to terms linear and quadratic in m s and linear in Am and m s Am reads: 



:) = jo- - r 2 Y - {a - r 2 )D 88 + -^=(n - r 2 ) J2 AuXi j m s , 
H Sl = I {r 2 K 2 (l - D 2 88 ) + (nKt - r 2 K 2 ) j ml 

ff dyn = ~\ i(No + 3iV 2 ) - 2N D 88 + (N - W 2 )D 2 S8 + 3(JVi - N 2 ) £ D 2 8 Jj 



h {1) = \-r 2 T 3 - ^(<r- r 2 )D 38 + (n - r 2 ) ^ ZW^l Am, 



h {2) 
"dyn 



^= |iVo£>38(3iV 2 - N )D 38 D 88 + 3{N 2 - JVi) ^ AjsAu j m s Am, 



CI = 3 | -r2K 2 D 38 D 8 s + {nK x - r 2 K 2 ) ^DasAw j m s Am, (14) 

where = Ki/Ii (i = 1,2) and T 3 stands for isospin. We have split the 0{m 2 ) and 0(m s Am) hamiltonian into the 
kinematical part which appears due to the fact that the quantization relation between the angular velocity and the 
spin operators J is m s and Am dependent (see Eq.(jll|)) and the dynamical part which comes from the expansion of 
the effective action in terms of m. 

The hamiltonians and ft/ 1 -* mix states of different SU(3) representations. The corresponding Oirn 2 ) contribution 
to the energy reads: 

' 5 -X-VV(3a + r 1 -4r 2 ) 2 l/ 2 m & 2 (15) 



250 V 2 2 



for the octet and for the decuplet: 



+ 356 (' " \ Y ~ l Y ') { ° + " ~ 2nf } hm '- (W) 

Here X = 1 — T(T + 1) + 1/4 Y 2 is the usual combination entering Gell-Mann-Okubo mass relations (T stands for 
isospin). 

Similarly one can write the 0(m s Am) wave function contribution to the octet states: 

e$ = i(<r - n) 2 (1 + Y) T 3 I 2 m s Am 

~^k {a + l ri ~ if 2)2 Y 713 /2TOsAm - 



4 



IV. HYPERON SPLITTINGS 



With the help of the matrix elements of the D functions and spin operators one arrives at the following result for 
the hyperon splittings: 

AM' 8 ' = A--Y-^=X-GY 2 , 
2 v/5 

AM (10) =B- —= Y-H Y 2 . (18) 
2y 2 

Constants A and B do not contribute to the splittings within the multiplets, however they shift the mass of the centers 
and contribute to the 10-8 mass difference. Constants G and H, not present in the first order Gell-Mann-Okubo 
mass formula, are of the order of m 2 . 

Hyperon splittings obtained with the help of Eq.(|l^) have been discussed in detail in Ref. Here for completeness 
we repeat only the main points. 

Experimentally one gets: 

F = E - N = 379 MeV, 
D = — (E — A) = 86 MeV, 

G=i(3A + S)-i(N + S)= 6.75 MeV (19) 

for the octet. For the decuplet the three operators: 1, Y and Y 2 do not form a complete basis and therefore there are 
two independent relations which determine constants C and H with small uncertainty: 

C = V2{E* - A) = -^=(ft - 2A + XT) = 422.5 ± 3.5 MeV, 
v2 

H = i(2£* - H* — A) = i(3£* - 2A - ft) = 2.83 ± 0.33 MeV. (20) 

In Tab. || we list the coefficients A ... H for a typical value of m s = 180 MeV as functions of the constituent mass 
M. It can be seen that in order to reproduce the experimental numbers of Eqs.(|l^,|2^) one has to take the constituent 
mass of the order of 400 MeV. Then all constants A. . . H are roughly reproduced. The constant G and H being 
of the order 0(m 2 ) are small. For reasonable strange quark masses 0(m 2 ) corrections to A,B,C and D are of the 
order of 20% of the leading 0(m s ) terms with the exception of F for which 0(m 2 ) corrections are almost zero. This 
is illustrated in Fig.l where 0(m 2 )/0(m s ) contributions to constants A...F are plotted as functions of m s for the 
fixed value of M = 420 MeV. 

TABLE II. Different contributions to coefficients of Eqs.(18) for M = 423 MeV and m s = 180 MeV 





0{m s ) 


kin. 


0{m 2 s ) 

dyn. 


w.f. 


total 


total 


exp. 


A 


546.10 


10.94 


-64.64 


-0.15 


-53.85 


492.25 




B 


546.10 


10.85 


-64.55 


-53.81 


-107.50 


438.60 




F 


381.20 


1.18 


-27.67 


22.76 


-3.73 


377.47 


379.00 


D 


120.76 


-0.02 


-11.78 


-0.07 


-11.87 


108.89 


86.00 


C 


348.16 


1.20 


-19.97 


90.92 


72.15 


420.32 


422.00 


G 


0.00 


0.61 


-0.66 


1.53 


1.48 


1.48 


6.75 


H 


0.00 


-0.29 


0.41 


4.57 


4.70 


4.70 


2.83 
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M=420 MeV 
0.3 i 



0.2 - 




50 100 150 200 250 

m s [MeV] 



FIG. 1. Ratios of 0(m^)/0(m B ) contributions to A ... D, F as functions of m s 
In order to make phenomcnological statements we adopt the following procedure: first for given M we find the 

—(8) 

optimal m s which reproduces 10-8 splitting. To this end we define the mean octet and decuplet values: M = 
1/2 (A + £) = 1155 MeV and M (10) = £* = 1385 MeV. Then A 10 _ 8 = M (1 ° } - M (8) = 230 MeV is given by: 

Aio-8 = 7^- + -B - A (21) 

All 

Since A — B = const, x ml one can numerically solve Eq. (pl|) for m s . The result is plotted in Fig. 2 (see also Tab. III). 



450 




390 1 

50 100 150 200 250 

m s [MeV] 



FIG. 2. M vs. m s dependence induced by the condition Aio-s = 230 MeV 
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In Fig. 3 we show the m s dependence of the deviations theory - experiment for each hyperon. One should remember 
that for each m s the optimal constituent quark mass M was used, so that Aio-s was automatically reproduced for 
each m s . The smallest deviations ±7 MeV for all splittings correspond to m s ~ 185 MeV, i.e. M ~ 426 MeV. 




FIG. 3. The m s dependence of the deviations theory - experiment for each hyperon for M chosen according to Fig. 2 

It is constructive to plot the dependence of constants F, D and C entering the Gell-Mann-Okubo mass relations 
as functions of m s (and optimal M). Fig. 4 shows the influence of the 0{m\) terms on these quantities. Dashed lines 
represent the first order results and solid lines the full result. It can be seen that 0(m%) corrections substantially 
improve model predictions. F hits its experimental value for m s = 181 MeV and the corresponding M ~ 424 MeV. 
For these values C is within the experimental errors and D overshoots slightly the experimental value. 

In this way we have fixed M - the only remaining free parameter of the model to the value of about 425 MeV. The 
corresponding value of m s ~ 180 MeV required to fit the 10-8 splitting overshoots the value deducted from the meson 
sector m s ~ 150 MeV. However, in an exact treatment of the vacuum sector and perturbation theory around this 
exact vacuum for the soliton sector, the crucial quantity is M s — M u , the difference of the constituent quark masses, 
instead of m s and this turns out to be ~ 180 — 200MeV ^> m s . Within such a treatment the discrepancy between 
meson and baryon sectors disappears E^ . 



TABLE III. Constants / and d as functions of M and m s chosen according to Fig. 2 



M 


m B 


f m 


J wf 


r(2) 
J dyn 


f (2) 
J kin 


£/ 




d (2) 

"wf 


d {2) 

"dyn 


d (2) 

"kin 


Ed 


395. 


19. 


3.34 


-0.05 


-0.01 


-0.00 


3.28 


0.46 


0.01 


-0.00 


-0.001 


0.47 


419. 


157. 


3.21 


-0.34 


-0.12 


-0.02 


2.72 


0.45 


0.07 


-0.02 


-0.006 


0.49 


423. 


177. 


3.18 


-0.37 


-0.14 


-0.02 


2.65 


0.45 


0.07 


-0.02 


-0.007 


0.49 


428. 


209. 


3.11 


-0.41 


-0.16 


-0.02 


2.53 


0.44 


0.08 


-0.03 


-0.007 


0.48 


442. 


258. 


3.03 


-0.45 


-0.19 


-0.02 


2.37 


0.43 


0.09 


-0.03 


-0.008 


0.48 
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optimal M 




150 
m s [MeV] 



250 



FIG. 4. F, C and D as functions of m s ; dashed lines - first order contributions, solid lines - full 0(m,g) result, note almost 
zero correction to F 



V. ISOSPIN SPLITTINGS 



Hadronic parts of the isospin splittings have been estimated in fig] : 

(n-p)n = 2.05 ±0.3, (XT - X+) H = 7.89 ± 0.3, (ET - S°) H = 5.5 ± 0.7 



(22) 



(in MeV) for the octet states. In Ref. || we have also estimated the splittings for the decuplet on the basis of the 
simple Dashen parametrization, which has proven to work equally well for both decuplet and octet particles. For the 
purpose of this work, however, we will concentrate entirely on the octet, since the numbers of Eq. (p2]) are experimental, 
whereas the decuplet estimations, in view of the lack of the data, were based upon the theoretical guess, as mentioned 
above. 

From theoretical point of view the isospin splittings are described by the formula analogous to the Gell-Mann-Okubo 
mass formula for the hyperon splittings, namely ff: 



(AM) H = -^fT 3 + dYT 3 , 

where, taking into account Eq. (|22|) , one gets: 

/ = 11.33 ± 1.14 MeV, d = 1.73 ± 0.38 MeV. 
On the other hand / and d can be directly evaluated in the present model with the help of Eqs.( fl4| , |l7| ): 



= \{? + n + 2r 2 ), = ^(2r 2 K 2 - 3n/f x ) m s 

(2) 



f (2) 



& = "i(2^o - m + 2N 2 ) m s , = --{a - r x f I 2 m s , 



1 



5 
S 1 



10 



3 

20 



(a-3r 1 + 2r 2 ), dg> = - 



(2r 2 K 2 - 5riJfi) m s 



7 (2) 

x dyn 



1 

45 



(AN - 5JV1 + 2N 2 ) m s , d. 



(2) 
wf 



1 / ^2 2 
— (a - ri) | a 

m y ' 125 



fx 4r 2 

y _ it 



h m s 



and 



/ = + /Si + /g, + OAm, d = (d« + C + C + «) A 



m. 



(23) 



(24) 



(25) 



(26) 
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It is worth noting that in the order m s Am the Gell-Mann-Okubo mass formula of Eq.(g3|) is not violated, but the 
coefficients / and d depend on m s . In Tab. Ill we list the values of / and d for different constituent masses M and 
for m s chosen according to Fig 2. 



It is clearly seen from Tab. [II that the second order corrections to / and d arc dominated by the wave function 
contribution, however the dynamical part, although smaller, is by no means negligible. This is in contrast to the 
hadronic splittings where the dynamical corrections were equally important as the wave function ones 

Finally in Figs. 5-7 we show the Am dependence of the isospin splittings for three constituent masses M = 419, 
423 and 428 MeV and three corresponding m s = 157, 177 and 209 MeV respectively. The first set corresponds to the 
strange quark mass as required by the meson sector. For this values, however, the hyperon spectrum is not correctly 
reproduced (see Figs. 3 and 4). The second set corresponds to the best fit to the hyperon spectra. The third set has 
been chosen to show what happens if m s starts overshooting the optimal value. 




M=423 MeV, m s =177 MeV 



-'V 




■ s - ' / 
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FIG. 6. Same as Fig.5 for M = 423 MeV and m s = 177 MeV 
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M=428 MeV, m s =209 MeV 
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m d -m u [MeV] 

FIG. 7. Same as Fig.5 for M = 428 MeV and m s = 209 MeV 

The figures clearly show, that the slopes of the splittings as functions of Am uniformly decrease with m s . Without 
the 0(m s ) corrections all splittings fall within the experimantal error bars for Am ~ 3.5 MeV (dashed lines on Figs. 
5-7). When the 0{m s ) corrections are included the Am range for which theoretical curves fall into experimental error 
bars shifts towards higher values and shrinks at the same time, so that for M = 428 MeV (Fig. 7) there is no common 
value of Am which would describe all splittings. For M — 423 MeV (Fig. 6) it is still possible to desribe all splittings 
wiith Am ~ 4.4 MeV. This value should be compared with the Am = 2.6 — 4.2 MeV which is needed to reproduce 
meson masses. 



VI. CONCLUDING REMARKS 

Almost all chiral models have problems in reproducing isospin splittings (see Ref. [p| and references therein). In 
our previous work || we have shown that in the semibosized NJL all isospin splitting within the octet and decuplet 
of baryons are reproduced within surprisingly good accuracy in the zeroth order in m s . In the Skyrme model the 
0(m s ) corrections to the isospin splittings are large and decrease the splittings. It has been therefore of importance 
to see how this corrections influence the mass differences in the present model. We have shown that the slope of the 
isospin splittings as functions of Am decreses with increasing m s . This is in principle what happens also in the SU(3) 
Skyrme model, however here the effect is much less pronounced. This is due to the existence of the anomalous parts 
associated with the quantities Ki$. Already at the zeroth order, where /W = 3.18 for M = 423 MeV the anomalous 



contribution amounts to 28%. The wave function corrections decrease / as seen from Tab. Ill by —0.37, however 
without the anomalous part the decrease would be much stronger: —0.46 for m s = 177 MeV. Altogether the final 
value of / = 2.65 would be decreased to 1.68 (i.e. by 37%) if the anomalous part was absent. On the contrary, the 
influence of anomalous terms on the coefficient d is small; it amounts to 4%. 

Let us stress once more that the successful phenomenology of the isospin splittings could not be achieved within the 
SU(2) NJL model, where, similarly to the Skyrme model, the matrix elements of the relevant operator simply vanish. 
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